A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network.

BACKGROUND: The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice. METHODS: FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription. RESULTS: FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency. CONCLUSIONS: Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

An Anterior Second Heart Field Enhancer Regulates the Gene Regulatory Network of the Cardiac Outflow Tract.

BACKGROUND: Conotruncal defects due to developmental abnormalities of the outflow tract (OFT) are an important cause of cyanotic congenital heart disease. Dysregulation of transcriptional programs tuned by NKX2-5 (NK2 homeobox 5), GATA6 (GATA binding protein 6), and TBX1 (T-box transcription factor 1) have been implicated in abnormal OFT morphogenesis. However, there remains no consensus on how these transcriptional programs function in a unified gene regulatory network within the OFT. METHODS: We generated mice harboring a 226-nucleotide deletion of a highly conserved cardiac enhancer containing 2 GATA-binding sites located ≈9.4 kb upstream of the transcription start site of Nkx2-5 (Nkx2-5∆enh) using CRISPR-Cas9 gene editing and assessed phenotypes. Cardiac defects in Nkx2-5∆enh/∆enh mice were structurally characterized using histology and scanning electron microscopy, and physiologically assessed using electrocardiography, echocardiography, and optical mapping. Transcriptome analyses were performed using RNA sequencing and single-cell RNA sequencing data sets. Endogenous GATA6 interaction with and activity on the NKX2-5 enhancer was studied using chromatin immunoprecipitation sequencing and transposase-accessible chromatin sequencing in human induced pluripotent stem cell-derived cardiomyocytes. RESULTS: Nkx2-5∆enh/∆enh mice recapitulated cyanotic conotruncal defects seen in patients with NKX2-5, GATA6, and TBX1 mutations. Nkx2-5∆enh/∆enh mice also exhibited defects in right Purkinje fiber network formation, resulting in right bundle-branch block. Enhancer deletion reduced embryonic Nkx2-5 expression selectively in the right ventricle and OFT of mutant hearts, indicating that enhancer activity is localized to the anterior second heart field. Transcriptional profiling of the mutant OFT revealed downregulation of important genes involved in OFT rotation and septation, such as Tbx1, Pitx2, and Sema3c. Endogenous GATA6 interacted with the highly conserved enhancer in human induced pluripotent stem cell-derived cardiomyocytes and in wild-type mouse hearts. We found critical dose dependency of cardiac enhancer accessibility on GATA6 gene dosage in human induced pluripotent stem cell-derived cardiomyocytes. CONCLUSIONS: Our results using human and mouse models reveal an essential gene regulatory network of the OFT that requires an anterior second heart field enhancer to link GATA6 with NKX2-5-dependent rotation and septation gene programs.

The cGAS-STING pathway is dispensable in a mouse model of LMNA-cardiomyopathy despite nuclear envelope rupture.

Mutations in the nuclear Lamin A/C gene (LMNA) cause diverse degenerative disorders, including malignant dilated cardiomyopathy in adults. A prevailing hypothesis postulates that LMNA mutations cause nuclear envelope ruptures that trigger pathogenic inflammatory signaling via the cGAS-STING cytosolic DNA-sensing pathway. Here, we provide evidence against this hypothesis, using a mouse model of LMNA-related cardiomyopathy that mimics Lamin A/C protein reduction observed in patient cardiomyocytes. We observed that pervasive nuclear envelope ruptures preceded the onset of cardiac transcriptional modulation and dilated cardiomyopathy. Nuclear ruptures activated DNA damage response without causing immediate cardiomyocyte death. However, cGAS-STING downstream cytokine genes remained inactive in the mutant cardiomyocytes. Deleting cGas or Sting did not alleviate cardiomyopathy. Instead, extracellular matrix signaling was predicted to emanate from Lamin A/C-reduced cardiomyocytes to communicate with fibroblasts in the heart. These findings suggest that cGAS-STING is not a major pathogenetic contributor to LMNA-related dilated cardiomyopathy in adult humans.

Single-cell genomics improves the discovery of risk variants and genes of atrial fibrillation.

Genome-wide association studies (GWAS) have linked hundreds of loci to cardiac diseases. However, in most loci the causal variants and their target genes remain unknown. We developed a combined experimental and analytical approach that integrates single cell epigenomics with GWAS to prioritize risk variants and genes. We profiled accessible chromatin in single cells obtained from human hearts and leveraged the data to study genetics of Atrial Fibrillation (AF), the most common cardiac arrhythmia. Enrichment analysis of AF risk variants using cell-type-resolved open chromatin regions (OCRs) implicated cardiomyocytes as the main mediator of AF risk. We then performed statistical fine-mapping, leveraging the information in OCRs, and identified putative causal variants in 122 AF-associated loci. Taking advantage of the fine-mapping results, our novel statistical procedure for gene discovery prioritized 46 high-confidence risk genes, highlighting transcription factors and signal transduction pathways important for heart development. In summary, our analysis provides a comprehensive map of AF risk variants and genes, and a general framework to integrate single-cell genomics with genetic studies of complex traits.

ETV2 primes hematoendothelial gene enhancers prior to hematoendothelial fate commitment.

Mechanisms underlying distinct specification, commitment, and differentiation phases of cell fate determination remain undefined due to difficulties capturing these processes. Here, we interrogate the activity of ETV2, a transcription factor necessary and sufficient for hematoendothelial differentiation, within isolated fate intermediates. We observe transcriptional upregulation of Etv2 and opening of ETV2-binding sites, indicating new ETV2 binding, in a common cardiac-hematoendothelial progenitor population. Accessible ETV2-binding sites are active at the Etv2 locus but not at other hematoendothelial regulator genes. Hematoendothelial commitment coincides with the activation of a small repertoire of previously accessible ETV2-binding sites at hematoendothelial regulators. Hematoendothelial differentiation accompanies activation of a large repertoire of new ETV2-binding sites and upregulation of hematopoietic and endothelial gene regulatory networks. This work distinguishes specification, commitment, and sublineage differentiation phases of ETV2-dependent transcription and suggests that the shift from ETV2 binding to ETV2-bound enhancer activation, not ETV2 binding to target enhancers, drives hematoendothelial fate commitment.

Tbx5 maintains atrial identity by regulating an atrial enhancer network.

Understanding how the atrial and ventricular chambers of the heart maintain their distinct identity is a prerequisite for treating chamber-specific diseases. Here, we selectively inactivated the transcription factor Tbx5 in the atrial working myocardium of the neonatal mouse heart to show that it is required to maintain atrial identity. Atrial Tbx5 inactivation downregulated highly chamber specific genes such as Myl7 and Nppa , and conversely, increased the expression of ventricular identity genes including Myl2 . Using combined single nucleus transcriptome and open chromatin profiling, we assessed genomic accessibility changes underlying the altered atrial identity expression program, identifying 1846 genomic loci with greater accessibility in control atrial cardiomyocytes compared to KO aCMs. 69% of the control-enriched ATAC regions were bound by TBX5, demonstrating a role for TBX5 in maintaining atrial genomic accessibility. These regions were associated with genes that had higher expression in control aCMs compared to KO aCMs, suggesting they act as TBX5-dependent enhancers. We tested this hypothesis by analyzing enhancer chromatin looping using HiChIP and found 510 chromatin loops that were sensitive to TBX5 dosage. Of the loops enriched in control aCMs, 73.7% contained anchors in control-enriched ATAC regions. Together, these data demonstrate a genomic role for TBX5 in maintaining the atrial gene expression program by binding to atrial enhancers and preserving tissue-specific chromatin architecture of atrial enhancers.

Frem1 activity is regulated by Sonic hedgehog signaling in the cranial neural crest mesenchyme during midfacial morphogenesis.

BACKGROUND: Frem1 has been linked to human face shape variation, dysmorphology, and malformation, but little is known about its regulation and biological role in facial development. RESULTS: During midfacial morphogenesis in mice, we observed Frem1 expression in the embryonic growth centers that form the median upper lip, nose, and palate. Expansive spatial gradients of Frem1 expression in the cranial neural crest cell (cNCC) mesenchyme of these tissues suggested transcriptional regulation by a secreted morphogen. Accordingly, Frem1 expression paralleled that of the conserved Sonic Hedgehog (Shh) target gene Gli1 in the cNCC mesenchyme. Suggesting direct transcriptional regulation by Shh signaling, we found that Frem1 expression is induced by SHH ligand stimulation or downstream pathway activation in cNCCs and observed GLI transcription factor binding at the Frem1 transcriptional start site during midfacial morphogenesis. Finally, we found that FREM1 is sufficient to induce cNCC proliferation in a concentration-dependent manner and that Shh pathway antagonism reduces Frem1 expression during pathogenesis of midfacial hypoplasia. CONCLUSIONS: By demonstrating that the Shh signaling pathway regulates Frem1 expression in cNCCs, these findings provide novel insight into the mechanisms underlying variation in midfacial morphogenesis.

Transcriptional Dysregulation Underlies Both Monogenic Arrhythmia Syndrome and Common Modifiers of Cardiac Repolarization.

BACKGROUND: Brugada syndrome (BrS) is an inherited arrhythmia syndrome caused by loss-of-function variants in the cardiac sodium channel gene SCN5A (sodium voltage-gated channel alpha subunit 5) in ≈20% of subjects. We identified a family with 4 individuals diagnosed with BrS harboring the rare G145R missense variant in the cardiac transcription factor TBX5 (T-box transcription factor 5) and no SCN5A variant. METHODS: We generated induced pluripotent stem cells (iPSCs) from 2 members of a family carrying TBX5-G145R and diagnosed with Brugada syndrome. After differentiation to iPSC-derived cardiomyocytes (iPSC-CMs), electrophysiologic characteristics were assessed by voltage- and current-clamp experiments (n=9 to 21 cells per group) and transcriptional differences by RNA sequencing (n=3 samples per group), and compared with iPSC-CMs in which G145R was corrected by CRISPR/Cas9 approaches. The role of platelet-derived growth factor (PDGF)/phosphoinositide 3-kinase (PI3K) pathway was elucidated by small molecule perturbation. The rate-corrected QT (QTc) interval association with serum PDGF was tested in the Framingham Heart Study cohort (n=1893 individuals). RESULTS: TBX5-G145R reduced transcriptional activity and caused multiple electrophysiologic abnormalities, including decreased peak and enhanced "late" cardiac sodium current (INa), which were entirely corrected by editing G145R to wild-type. Transcriptional profiling and functional assays in genome-unedited and -edited iPSC-CMs showed direct SCN5A down-regulation caused decreased peak INa, and that reduced PDGF receptor (PDGFRA [platelet-derived growth factor receptor α]) expression and blunted signal transduction to PI3K was implicated in enhanced late INa. Tbx5 regulation of the PDGF axis increased arrhythmia risk due to disruption of PDGF signaling and was conserved in murine model systems. PDGF receptor blockade markedly prolonged normal iPSC-CM action potentials and plasma levels of PDGF in the Framingham Heart Study were inversely correlated with the QTc interval (P<0.001). CONCLUSIONS: These results not only establish decreased SCN5A transcription by the TBX5 variant as a cause of BrS, but also reveal a new general transcriptional mechanism of arrhythmogenesis of enhanced late sodium current caused by reduced PDGF receptor-mediated PI3K signaling.

Tbx5 maintains atrial identity in post-natal cardiomyocytes by regulating an atrial-specific enhancer network.

Understanding how the atrial and ventricular heart chambers maintain distinct identities is a prerequisite for treating chamber-specific diseases. Here, we selectively knocked out (KO) the transcription factor Tbx5 in the atrial working myocardium to evaluate its requirement for atrial identity. Atrial Tbx5 inactivation downregulated atrial cardiomyocyte (aCM) selective gene expression. Using concurrent single nucleus transcriptome and open chromatin profiling, genomic accessibility differences were identified between control and Tbx5 KO aCMs, revealing that 69% of the control-enriched ATAC regions were bound by TBX5. Genes associated with these regions were downregulated in KO aCMs, suggesting they function as TBX5-dependent enhancers. Comparing enhancer chromatin looping using H3K27ac HiChIP identified 510 chromatin loops sensitive to TBX5 dosage, and 74.8% of control-enriched loops contained anchors in control-enriched ATAC regions. Together, these data demonstrate TBX5 maintains the atrial gene expression program by binding to and preserving the tissue-specific chromatin architecture of atrial enhancers.

Hedgehog signaling activates a mammalian heterochronic gene regulatory network controlling differentiation timing across lineages.

Many developmental signaling pathways have been implicated in lineage-specific differentiation; however, mechanisms that explicitly control differentiation timing remain poorly defined in mammals. We report that murine Hedgehog signaling is a heterochronic pathway that determines the timing of progenitor differentiation. Hedgehog activity was necessary to prevent premature differentiation of second heart field (SHF) cardiac progenitors in mouse embryos, and the Hedgehog transcription factor GLI1 was sufficient to delay differentiation of cardiac progenitors in vitro. GLI1 directly activated a de novo progenitor-specific network in vitro, akin to that of SHF progenitors in vivo, which prevented the onset of the cardiac differentiation program. A Hedgehog signaling-dependent active-to-repressive GLI transition functioned as a differentiation timer, restricting the progenitor network to the SHF. GLI1 expression was associated with progenitor status across germ layers, and it delayed the differentiation of neural progenitors in vitro, suggesting a broad role for Hedgehog signaling as a heterochronic pathway.

Differential Etv2 threshold requirement for endothelial and erythropoietic development.

Endothelial and erythropoietic lineages arise from a common developmental progenitor. Etv2 is a master transcriptional regulator required for the development of both lineages. However, the mechanisms through which Etv2 initiates the gene-regulatory networks (GRNs) for endothelial and erythropoietic specification and how the two GRNs diverge downstream of Etv2 remain incompletely understood. Here, by analyzing a hypomorphic Etv2 mutant, we demonstrate different threshold requirements for initiation of the downstream GRNs for endothelial and erythropoietic development. We show that Etv2 functions directly in a coherent feedforward transcriptional network for vascular endothelial development, and a low level of Etv2 expression is sufficient to induce and sustain the endothelial GRN. In contrast, Etv2 induces the erythropoietic GRN indirectly via activation of Tal1, which requires a significantly higher threshold of Etv2 to initiate and sustain erythropoietic development. These results provide important mechanistic insight into the divergence of the endothelial and erythropoietic lineages.

Detecting critical transition signals from single-cell transcriptomes to infer lineage-determining transcription factors.

Analyzing single-cell transcriptomes promises to decipher the plasticity, heterogeneity, and rapid switches in developmental cellular state transitions. Such analyses require the identification of gene markers for semi-stable transition states. However, there are nontrivial challenges such as unexplainable stochasticity, variable population sizes, and alternative trajectory constructions. By advancing current tipping-point theory-based models with feature selection, network decomposition, accurate estimation of correlations, and optimization, we developed BioTIP to overcome these challenges. BioTIP identifies a small group of genes, called critical transition signal (CTS), to characterize regulated stochasticity during semi-stable transitions. Although methods rooted in different theories converged at the same transition events in two benchmark datasets, BioTIP is unique in inferring lineage-determining transcription factors governing critical transition. Applying BioTIP to mouse gastrulation data, we identify multiple CTSs from one dataset and validated their significance in another independent dataset. We detect the established regulator Etv2 whose expression change drives the haemato-endothelial bifurcation, and its targets together in CTS across three datasets. After comparing to three current methods using six datasets, we show that BioTIP is accurate, user-friendly, independent of pseudo-temporal trajectory, and captures significantly interconnected and reproducible CTSs. We expect BioTIP to provide great insight into dynamic regulations of lineage-determining factors.

Tbx5 drives Aldh1a2 expression to regulate a RA-Hedgehog-Wnt gene regulatory network coordinating cardiopulmonary development.

The gene regulatory networks that coordinate the development of the cardiac and pulmonary systems are essential for terrestrial life but poorly understood. The T-box transcription factor Tbx5 is critical for both pulmonary specification and heart development, but how these activities are mechanistically integrated remains unclear. Here using Xenopus and mouse embryos, we establish molecular links between Tbx5 and retinoic acid (RA) signaling in the mesoderm and between RA signaling and sonic hedgehog expression in the endoderm to unveil a conserved RA-Hedgehog-Wnt signaling cascade coordinating cardiopulmonary (CP) development. We demonstrate that Tbx5 directly maintains expression of aldh1a2, the RA-synthesizing enzyme, in the foregut lateral plate mesoderm via an evolutionarily conserved intronic enhancer. Tbx5 promotes posterior second heart field identity in a positive feedback loop with RA, antagonizing a Fgf8-Cyp regulatory module to restrict FGF activity to the anterior. We find that Tbx5/Aldh1a2-dependent RA signaling directly activates shh transcription in the adjacent foregut endoderm through a conserved MACS1 enhancer. Hedgehog signaling coordinates with Tbx5 in the mesoderm to activate expression of wnt2/2b, which induces pulmonary fate in the foregut endoderm. These results provide mechanistic insight into the interrelationship between heart and lung development informing CP evolution and birth defects.

Identification of direct transcriptional targets of NFATC2 that promote ß cell proliferation.

The transcription factor NFATC2 induces ß cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets. By integrating changes in gene expression with genomic binding sites for NFATC2, we identified approximately 2200 transcriptional targets of NFATC2. Genes induced by NFATC2 were enriched for transcripts that regulate the cell cycle and for DNA motifs associated with the transcription factor FOXP. Islets from an endocrine-specific Foxp1, Foxp2, and Foxp4 triple-knockout mouse were less responsive to NFATC2-induced ß cell proliferation, suggesting the FOXP family works to regulate ß cell proliferation in concert with NFATC2. NFATC2 induced ß cell proliferation in both mouse and human islets, whereas NFATC1 did so only in human islets. Exploiting this species difference, we identified approximately 250 direct transcriptional targets of NFAT in human islets. This gene set enriches for cell cycle-associated transcripts and includes Nr4a1. Deletion of Nr4a1 reduced the capacity of NFATC2 to induce ß cell proliferation, suggesting that much of the effect of NFATC2 occurs through its induction of Nr4a1. Integration of noncoding RNA expression, chromatin accessibility, and NFATC2 binding sites enabled us to identify NFATC2-dependent enhancer loci that mediate ß cell proliferation.

Control of cardiomyocyte differentiation timing by intercellular signaling pathways.

Congenital Heart Disease (CHD), malformations of the heart present at birth, is the most common class of life-threatening birth defect (Hoffman (1995) [1], Gelb (2004) [2], Gelb (2014) [3]). A major research challenge is to elucidate the genetic determinants of CHD and mechanistically link CHD ontogeny to a molecular understanding of heart development. Although the embryonic origins of CHD are unclear in most cases, dysregulation of cardiovascular lineage specification, patterning, proliferation, migration or differentiation have been described (Olson (2004) [4], Olson (2006) [5], Srivastava (2006) [6], Dunwoodie (2007) [7], Bruneau (2008) [8]). Cardiac differentiation is the process whereby cells become progressively more dedicated in a trajectory through the cardiac lineage towards mature cardiomyocytes. Defects in cardiac differentiation have been linked to CHD, although how the complex control of cardiac differentiation prevents CHD is just beginning to be understood. The stages of cardiac differentiation are highly stereotyped and have been well-characterized (Kattman et al. (2011) [9], Wamstad et al. (2012) [10], Luna-Zurita et al. (2016) [11], Loh et al. (2016) [12], DeLaughter et al. (2016) [13]); however, the developmental and molecular mechanisms that promote or delay the transition of a cell through these stages have not been as deeply investigated. Tight temporal control of progenitor differentiation is critically important for normal organ size, spatial organization, and cellular physiology and homeostasis of all organ systems (Raff et al. (1985) [14], Amthor et al. (1998) [15], Kopan et al. (2014) [16]). This review will focus on the action of signaling pathways in the control of cardiomyocyte differentiation timing. Numerous signaling pathways, including the Wnt, Fibroblast Growth Factor, Hedgehog, Bone Morphogenetic Protein, Insulin-like Growth Factor, Thyroid Hormone and Hippo pathways, have all been implicated in promoting or inhibiting transitions along the cardiac differentiation trajectory. Gaining a deeper understanding of the mechanisms controlling cardiac differentiation timing promises to yield insights into the etiology of CHD and to inform approaches to restore function to damaged hearts.

Identifying genetic factors that contribute to the increased risk of congenital heart defects in infants with Down syndrome.

Atrioventricular septal defects (AVSD) are a severe congenital heart defect present in individuals with Down syndrome (DS) at a > 2000-fold increased prevalence compared to the general population. This study aimed to identify risk-associated genes and pathways and to examine a potential polygenic contribution to AVSD in DS. We analyzed a total cohort of 702 individuals with DS with or without AVSD, with genomic data from whole exome sequencing, whole genome sequencing, and/or array-based imputation. We utilized sequence kernel association testing and polygenic risk score (PRS) methods to examine rare and common variants. Our findings suggest that the Notch pathway, particularly NOTCH4, as well as genes involved in the ciliome including CEP290 may play a role in AVSD in DS. These pathways have also been implicated in DS-associated AVSD in prior studies. A polygenic component for AVSD in DS has not been examined previously. Using weights based on the largest genome-wide association study of congenital heart defects available (2594 cases and 5159 controls; all general population samples), we found PRS to be associated with AVSD with odds ratios ranging from 1.2 to 1.3 per standard deviation increase in PRS and corresponding liability r2 values of approximately 1%, suggesting at least a small polygenic contribution to DS-associated AVSD. Future studies with larger sample sizes will improve identification and quantification of genetic contributions to AVSD in DS.

Chamber-specific transcriptional responses in atrial fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia, yet the molecular signature of the vulnerable atrial substrate is not well understood. Here, we delineated a distinct transcriptional signature in right versus left atrial cardiomyocytes (CMs) at baseline and identified chamber-specific gene expression changes in patients with a history of AF in the setting of end-stage heart failure (AF+HF) that are not present in heart failure alone (HF). We observed that human left atrial (LA) CMs exhibited Notch pathway activation and increased ploidy in AF+HF but not in HF alone. Transient activation of Notch signaling within adult CMs in a murine genetic model is sufficient to increase ploidy in both atrial chambers. Notch activation within LA CMs generated a transcriptomic fingerprint resembling AF, with dysregulation of transcription factor and ion channel genes, including Pitx2, Tbx5, Kcnh2, Kcnq1, and Kcnip2. Notch activation also produced distinct cellular electrophysiologic responses in LA versus right atrial CMs, prolonging the action potential duration (APD) without altering the upstroke velocity in the left atrium and reducing the maximal upstroke velocity without altering the APD in the right atrium. Our results support a shared human/murine model of increased Notch pathway activity predisposing to AF.

Epigenetic and Transcriptional Networks Underlying Atrial Fibrillation.

Genome-wide association studies have uncovered over a 100 genetic loci associated with atrial fibrillation (AF), the most common arrhythmia. Many of the top AF-associated loci harbor key cardiac transcription factors, including PITX2, TBX5, PRRX1, and ZFHX3. Moreover, the vast majority of the AF-associated variants lie within noncoding regions of the genome where causal variants affect gene expression by altering the activity of transcription factors and the epigenetic state of chromatin. In this review, we discuss a transcriptional regulatory network model for AF defined by effector genes in Genome-wide association studies loci. We describe the current state of the field regarding the identification and function of AF-relevant gene regulatory networks, including variant regulatory elements, dose-sensitive transcription factor functionality, target genes, and epigenetic states. We illustrate how altered transcriptional networks may impact cardiomyocyte function and ionic currents that impact AF risk. Last, we identify the need for improved tools to identify and functionally test transcriptional components to define the links between genetic variation, epigenetic gene regulation, and atrial function.

Hedgehog-FGF signaling axis patterns anterior mesoderm during gastrulation.

The mechanisms used by embryos to pattern tissues across their axes has fascinated developmental biologists since the founding of embryology. Here, using single-cell technology, we interrogate complex patterning defects and define a Hedgehog (Hh)-fibroblast growth factor (FGF) signaling axis required for anterior mesoderm lineage development during gastrulation. Single-cell transcriptome analysis of Hh-deficient mesoderm revealed selective deficits in anterior mesoderm populations, culminating in defects to anterior embryonic structures, including the pharyngeal arches, heart, and anterior somites. Transcriptional profiling of Hh-deficient mesoderm during gastrulation revealed disruptions to both transcriptional patterning of the mesoderm and FGF signaling for mesoderm migration. Mesoderm-specific Fgf4/Fgf8 double-mutants recapitulated anterior mesoderm defects and Hh-dependent GLI transcription factors modulated enhancers at FGF gene loci. Cellular migration defects during gastrulation induced by Hh pathway antagonism were mitigated by the addition of FGF4 protein. These findings implicate a multicomponent signaling hierarchy activated by Hh ligands from the embryonic node and executed by FGF signals in nascent mesoderm to control anterior mesoderm patterning.